Lattice dielectric response of CdCu3Ti40i2 and of CaCu3Ti40i2 from first principles 
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Structural, vibrational, and lattice dielectric properties of CdCu3Ti40i2 are studied using density- 
functional theory within the local spin-density approximation, and the results are compared with 
those computed previously for CaCu3Ti40i2. Replacing Ca with Cd is found to leave many calcu- 
lated quantities largely unaltered, although significant difi^erences do emerge in zone-center optical 
phonon frequencies and mode efi'ective charges. The computed phonon frequencies of CdCu3Ti40i2 
are found to be in excellent agreement with experiment, and the computed lattice contribution to the 
intrinsic static dielectric constant (~60) also agrees exceptionally well with a recent optical absorp- 
tion experiment. These results provide further support for a picture in which the lattice dielectric 
response is essentially conventional, suggesting an extrinsic origin for the anomalous low-frequency 
dielectric response recently observed in both materials. 



PACS: 77.84.-S, 63.20.-e, 71.20.-b, 77.22.-d 

The complex perovskite compound CaCu3Ti40i2 
(CCTO) has recently attracted considerable attention 
in view of its anomalously large dielectric response. 
Low-frequency studies rugcovered dielectric constants 
near 161)00 for ceramicsoa and 80,000 for single-crystal 
samplesH that are nearly constant over a wide tempera- 
ture range (~100-600K). Moreover, above a critical fre- 
quency that ranges between '--^lOHz and 1 MHz depend- 
ing on temperature, the dielectric constant drops to a 
rather ordinary value, of '^100, typical of conventional 
perovskite insulatoraj where contributions from zone- 
center infrared-active phonons can easily result in values 
of this magnitude. Detailed measurements of the dy- 
namic susceptibility indicated that the low-temperature 
crossover is characteristic of an activated Debye-type re- 
laxation process .0 p. 

In a previous first-principles study of CCTOja we cal- 
culated the lattice contribution to the dielectric constant, 
and found it to be in order-of-m£iffnitude agreement with 
the measured far-infrared valueu and drastically lower 
than the static value. Thus, while the dielectric response 
in the infrared range is essentially well understood, the 
origin of the enormous low-frequency response and its 
dynamic behavior remain perplexing. An extrinsic mech- 
anism, i.e., one associated with defects, domain bound- 
aries, or other crystalline deficiencies, was initially pro- 
posed by Subramanian et aln These authors reported 
twinning in CCTO samples and suggested that internal 
capacitive barrier layers, formed by domain boundaries, 
might explain the unusual dielectric response. A siiru 
ilar explanation has been proposed by another group.El 
We discussed some of these issues in Ref. || and, in the 
absence of direct evidence for intrinsic lattice or elec- 
tronic origins, we concluded that extrinsic effects were 
the likely source of the unusual low-frequencK behavior. 
A subsequent report by Lunkenheimer et al.u discussed 
a scenario in which barrier layers associated with poor 
metallic contacts are responsible for the Debye relax- 



ation. Further discjiission of extrinsic mechanisms will 
be given elsewhere. cl 

While the enormous low-frequency response is the pri- 
mary puzzle in this material, two other aspects of the 
response at higher frequencies are also currently enig- 
matic. First, the oscillator strength of one low fre- 
quency IR-active mode has been observed to increase 
strongly with decreasing temperature. And second, one 
IR-active mode, at somewhat higher frequency, is not 
observed at all, although it is expected by symmetry 
and has been predicted to possess a non-vanishing os- 
cillator strcngth.ErLi It is not clear whether there could 
be any connection between these IR anomalies and the 
enormous low-frequency response, but in any case, fur- 
ther experimental and theoretical scrutiny of the infrared 
phonon spectrum would clearly be of significant physical 
interest .d 

In this report, we extend our previous work on CCTO 
by carrying out a parallel study oii-the closely related ma- 
terial CdCu3Ti40i2 (CdCTO).^! Aj-ecent experimen- 
tal study of ceramic CdCTO samplesQ has revealed that 
CdCTO also has a temperature-independent anomalous 
response over a wide range and shows a similar De- 
bye crossover between low- and high-frequency regimes. 
While the magnitude of this response (e ~400 for ceramic 
samples) is considerable, it is nonetheless significantly 
lower than in CCTO, which is somewhat surprising since 
Cd and Ca are chemically so similar. (Anomalies also 
appear in the low-frequency IR-active mode oscillator 
strengths of CdCTO.) After calculating ground-state and 
dielectric properties of CdCTO, we compare with recent 
experiments and with the corresponding quantities for 
CCTO. The results obtained for both materials are sim- 
ilar, typical of ordinary perovskite oxides, and iiijexccl- 
lent agreement with recent optical experiments. mj Our 
results indicate no significant fundamental differences in 
the intrinsic properties of these two materials and thus 
provide support for an extrinsic origin of their substantial 
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TABLE I. Comparison of calculated and measured struc- 
tural parameters of CCTO and CdCTO. Both have the space 
group 7m3 (point group Th); the Wyckoff positions are 
Ca/Cd(0, 0, 0), Cu(l/2, 0, 0), Ti(l/4, 1/4, 1/4), 0(a;, y, 0). See 
Ref. 13 for further details. 





Structural parameter 


LSDA 


Exp 




a (A) 


7.324 


7.384^ 


CaCu3Ti40i2 


Oix) 


0.303 


0.303 




0{y) 


0.175 


0.179 




a (A) 


7.324 


7.384" 


CdCu3Ti40i2 


Oix) 


0.306 






0{y) 


0.176 





^35 K (a=7.39lA at 298 K). 
•"298 K. 



low-frequency responses. 

Our first-principles study of CdCTO employs density- 
functional theory within the local spin-density ap- 
proximation (LSDA), as implemented withiH|— the Vi- 
enna ab-initio Simulations Package (VASP).t3y VASP 
utilizes a plane-|W|ave basis and Vanderbilt ultrasoft 
pseudopotentialsll3 AlL_pseudopotentials include non- 
linear core corrections.lij For Ca and Ti, these poten- 
tials treat the highest occupied p shell electrons explic- 
itly as valence; for Cu, we also consider electrons in the 
3d and 4s shells self-consistently; and for Cd, the 4d and 
5s shells are included explicitly. The ions are relaxed 
towards equilibrium until the Hellmann-Feynman forces 
are less than 10^'^eV/A. Brillouin-zone integrations are 
performed with a Gaussian broadening of 0.1 eV during 
all relaxations. A 37 Ry plane-wave cut-off and a 2x2x2 
Monkhorst-Pack k-point mesh (equivalent to a 4x4x4 
mesh for a single 5-atom perovskite unit cell) results in 
good convergence of all properties resorted here. 

Neutron-diffraction measurementsEJ suggest that, as 
with CCTO, CdCTO crystallizes in a lattice with a 20- 
atom body-centered cubic primitive celllj having space 
and point groups /m3 and Th, respectively. Since the 
/m3 space group is centrosymmetric, the spontaneous 
polarization necessarily vanishes by symmetry. To ac- 
commodate the observed antiferromagnetic (AFM) spin 
arrangement, all calculations are performed using a dou- 
bled 40-atom simple-cubic unit cell containing eight per- 
ovskite 5-atom units, of which 3/4 have Cu on the A site 
and 1/4 have Ca or Cd on the A site. (Each Cu-Cu 
nearest- neighbor pair has antiparallel spins.) 

Structural parameters resulting from the relaxation ap- 
pear in Table I. Replacing Ca_with Cd leaves the lattice 
constant virtually unchangedjl3 and the internal param- 
eters also remain essentially the same. Only minor struc- 
tural changes would be expected, given the similar nom- 
inal valence (+2) aiid ionic radii of Ca and Cd (1.48 and 
1.45 A, respectivel5it3) . Unfortunately we are unaware of 
any experimental refinements of the internal structural 
parameters of CdCTO, and thus direct comparison of 
our computed oxygen positions with experiment is not 



possible. As for the lattice parameters, however, our 
LSDA results are consistent with the near-negligible dif- 
ference (<0.1%) between the lattice constants measured 
for CCTO and CdCTO (7.391 and 7.384 A at 298 K, 
respectively). (The CCTO lattice constant at 35 K is 
smaller and, coincidentally, 7.384 A; a low-temperature 
lattice constant for CdCTO has yet to be measured.) 
Lattice constants computed within the generalized gra- 
dient approximationim (GGA) are calculated to be about 
1.2% larger than experiment for both CCTO and Cd- 
CTO. All results presented below were computed within 
the LSDA without gradient corrections. 

The computed electronic structure (at T = 0) is quite 
similar for both materials; replacing Ca with Cd evi- 
dently has little effect on electronic and magnetic prop- 
erties. CdCTO is found to possess an antiferromagnetic 
(AFM) insulating ground state, and its single-particle 
density of states (DOS) closely corresponds to that of 
CCTO near the band edges (see Figs. 2 and 3 in Ref. ||). 
Although the Cd ion has a nominal valence of +2, its 
electronic configuration includes a filled Ad shell, which 
results in a well-localized set of d bands located about 
8 eV below the valence band maximum. As with CCTO, 
the gap and magnetic moment of CdCTO originate from 
a splitting of hybridized Cu(3d)-0(2p) cr-antibonding 
states. The magnetic moment of each CUO4 plaque- 
tte, estimated from spin densities, is ~0.84/iB, and the 
(indirect) band gap is computed to be about 0.19 eV. 
To our knowledge, experimental values for the optical 
gap remain unavailable. However, the observed gap will 
undoubtedly be strongly underestimated by our calcula- 
tions, as was found for CCTO, owing to the well-known 
limitations of the LSDA. 

The lattice contribution to the dielectric constant that 
we computed for CCTO ('^40) was roughly a factor 
of two smaller than the value measured by far-infrared 
spectroscopyJj We now perform a similar calculation for 
CdCTO and carefully examine any differences between 
the materials. 

The static lattice dielectric response can be approxi- 
mated (neglecting anharmonicity) as the zero-frequency 
response of a system of classical Lorentz oscillators, i.e., 

where uj\ and Z'^ are respectively the IR-active mode 
frequencies and mode dynamical charges, and flo^ = 
Aire'^ /mcjV is a characteristic frequency having the in- 
terpretation of a plasma frequency of a gas of objects of 
mass mo=lamu, charge e, and density (V is the 

20- atom primitive cell volume). 

The lattice dielectric constant Cph is obtained by first 
calculating the frequency and mode effective charges of 
each zone-center IR-active phonon, and then inserting 
these quantities into Eq. (1). Since CdCTO and CCTO 
possess the same structure, their zone-center phonons are 
computed and analyzed as in our previous studyU As 
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TABLE II. Comparison of calculated mode frequencies ljx, 
effective charges Z^, and oscillator strengths S\ = Z"^ Iu\ 
of IR-active T„ modes with the experimental values (at 
T = 295 K) of Homes et al. (Refs. |, §). See Ref. ^ for 
further details. 



CaCu3Ti40i2 CdCu3Ti40i2 



uj (cm 


Sa 




UJ (cm" 




Sa 




LSDA 


Exper. 


LSDA Exper. 


LSDA Exper. 


LSDA Exper. 


125 


122.3 


4.6 


14.3 


72 


74 


15.1 


22.8 


135 


140.8 


9.3 


15.9 


125 


121 


7.3 


11.3 


158 


160.8 


7.5 


6.92 


148 


155 


0.4 


7.9 


199 


198.9 


2.7 


5.25 


170 


167 


11.8 


4.3 


261 


253.9 


12.5 


13.8 


238 


239 


12.6 


7.9 


310 


307.6 


1.0 


0.68 


303 


295 


2.4 


2.8 


385 


382.1 


0.5 


1.96 


385 


385 


0.7 


1.6 


416 


421.0 


5.5 


1.72 


405 


422 


2.4 


1.4 


471 




1.5 




461 


469 


4.2 


0.3 


494 


504.2 


0.6 


0.78 


499 


494 


0.4 


1.1 


547 


552.4 


0.4 


0.62 


545 


550 


0.3 


0.7 



discussed there, of the six irreducible representations al- 
lowed by the Th point group, only the T„ modes display 
IR activity. After obtaining symmetry-adapted modes 
of T„ symmetry, IR-active phonon frequencies and cor- 
responding eigenvectors are calculated using the fcpzen- 
phonon method, as described in detail elsewhereE The 
results appear in Table II, where we also compare with 
corresponding theoretical results for CCTOEa and with 
recent experimental data. As was the case for CCTO, all 
CdCTO IR-active modes are found to be stable (w^ > 0) 
and their computed frequencies agree very well with 
experiments 

Given that Ca and Cd are isoelectronic, and that 
CCTO and CdCTO possess similar structural and elec- 
tronic properties, we may expect their phonon frequen- 
cies and eigenvectors to be similar as well, especially for 
modes that do not involve Ca or Cd. To investigate this 
and related issues, we find it useful to sort the zone- 
center mode frequencies from lowest to highest in fre- 
quency, and then project each CCTO eigenvector onto 
the corresponding one of CdCTO. For all but the lowest 
four modes, the projections result in values > 0.9, re- 
flecting a close correspondence between the character of 
their eigenmodes. However, the third- and fourth-lowest 
eigenvectors (the 148 and 170 cm~^ modes of CdCTO 
and the 158 and 199 cm~^ modes of CCTO) are substan- 
tially mixed by replacement of Ca with Cd, and possess 
considerably different eigenvectors. Further, the frequen- 
cies of the lowest two modes reverse order on going to the 
Cd from the Ca material: the 72cm~^ mode of CdCTO 
shares nearly the same eigenvector with the 141 cm^^ 
mode in CCTO, and we observe a close correspondence 
between the 125 cm"! modes of CdCTO and CCTO. The 
latter two are dominated by Cu-0 motion, and thus are 
unaffected by the replacement of Ca with Cd. We note 
here that peculiar changes in oscillator strengths with 



TABLE III. Comparison of calculated Born effective 
charge tensors Z* (in a Cartesian basis) for CCTO and Cd- 
CTO. Atoms are at Wyckoff positions given in the caption of 
Table I. 





CaCu3Ti40i2 




CdCu3Ti40l2 




Ua/Uu 


2.46 








2.47 





U 







2.46 








2.47 













2.46 





U 


2.47 


Cu 


2.06 








2.09 













1.85 








1.87 













1.21 








1.10 


Ti 


6.98 


-0.13 


0.06 


6.91 


-0.18 


0.03 




0.06 


6.98 - 


-0.13 


0.03 


6.91 - 


-0.18 




-0.13 


0.06 


6.98 


-0.18 


0.03 


6.91 





-1.92 


0.54 





-1.94 


0.54 







0.22 


-1.94 





0.25 


-1.91 













-5.01 








-4.94 



temperature are observed to occur for this mode in both 
materials, indirect confirmation of their similar charac- 
ter. 

The conspicuous reduction in frequency of the 
135cm-i mode in CCTO to 72cm-i in CdCTO origi- 
nates from both inertial and chemical effects. The masses 
of Ca and Cd are 112.4 and 40.1 amu respectively; on this 
basis alone, the frequency could decline at most by a ratio 
of y\Mca/A/cd =0.597. Since Ca/Cu motion accounts for 
only 38% of the eigenvector character, the inertial shift 
must actually be considerably smaller, and thus cannot 
account for the calculated frequencies. Indeed, we car- 
ried out a test calculation in which the Ca mass was 
replaced by that of Cd while retaining the force-constant 
matrix of CCTO. Rediagonalizing the dynamical matrix, 
we observe that the lowest four frequencies change to 102, 
126, 155, and 186 cm~^, respectively (the other seven re- 
main almost unchanged). Thus, the mass effect alone 
would only be able to explain a decrease from 135 cm~^ 
to 102cm^^, i.e., a downward shift of about 24%. This 
is much less than the reduction observed experimentally 
(48%) or theoretically (47%). The remainder of the shift 
is associated with an appreciable softening (~ 40%) of 
the force constants of this mode. Evidently this soften- 
ing is connected with chemical differences between the 
two atoms, e.g., the Ad shell present in Cd but absent in 

In our previous work,i we calculated the mode effective 
charges = VAPz/uq for each z-polarized lEsactive 
mode of CCTO using the Berry-phase formalism.llj (Here 
AP^ = Pz{ueq - uqUx) - Pz{ueq), Ug, Is the equilibrium 
configuration of ions, and ua is the eigenvector of mode 
A and uq is its amplitude.) We repeat the process for 
CdCTO and obtain the mode effective charges; the cor- 
responding oscillator strengths are compared with exper- 
iment in Table II. As we found for CCTO, the individual 
oscillator strengths disagree with experimental values for 
some of the modes. Remarkably, however, the mode that 
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was absent from the CCTO spectrum reappears in the 
case of the CdCTO spectrum (Ref. ||), albeit with an os- 
cillator strength of only 0.3; for both materials this mode 
is primarily composed of O 80%) and Ti (~ 15%) dis- 
placements. 

From the mode effective charges we now estimate 
the lattice contribution to the static lattice dielec- 
tric constant by summing over our calculated oscillator 
strengths, i.e., eph = Ea ^aV^a = Ea'S'a- Using 
the eigenmodes and their effective charges in Table II, 
we obtain Cph =57.6, in excellent agreement with the ex- 
perimental value at room temperature (cph =62.1) but 
somewhat less that the low temperature value (eph =106 
at 10 K), where an anomalous increase in the oscilla- 
tor strength has been observedE Compared with value 
of -45 computed for CCTO, that obtained for CdCTO 
is —50% larger, a consequence of the lowest-frequency IR 
phonon in CdCTO having a lower frequency and larger 
effective charge than the corresponding mode in CCTO. 
Thus, despite differences in computed oscillator strengths 
for these modes, their combined contribution to the di- 
electric constant is remarkably close to that observed ex- 
perimentally. Since the computed frequencies agree much 
better with experiment than the oscillator strengths, the 
agreement for eph implies some degree of error cancella- 
tion in the latter. 

To complete our comparative study of these two ma- 
terials, we also compute their atomic effective charge 
tensors. This is done in practice by starting with the 
mode effective charges and then using the correspond- 
ing mode eigenvectors to transform back to an atomic 
displacement basis. (The results for CCTO were not re- 
ported in Ref. || and are given for the first time here.) As 
can be seen from Table III, the atomic effective charge 
tensors are virtually identical for both materials,c2l and 
thus the same discussion applies to both. The Ca and 
Cd cations have diagonal and isotropic effective charge 
tensors of almost identical magnitude, 4-2.46 and -t-2.47 
respectively. Thus, while the force constants are quite 
sensitive to the different underlying atomic shell struc- 
ture of the Ca and Cd ions, the effective charge is not, 
suggesting that there is little dynamical charge transfer 
between this cation and its oxygen neighbors. For Cu, 
the two components associated with displacements par- 
allel to the Cu04 plaquettes are close to the nominal va- 
lence (-|-2), while the third component corresponding to 
normal displacements is considerably smaller. Except for 
the presence of some small off-diagonal components per- 
mitted by symmetry, the results for Ti and O are strongly 
reminiscent of other perovskites: we find that Z* for Ti 
has a large positive anomalous component (relative to 
the nominal -f 4) , and correspondingly each oxygen has a 
negative anomalous component (relative to the nominal 
—2) for displacements toward the nearest-neighbor Ti. 
Indeed, the atomic effective charges computed for Ca, 
Cd, Ti, and O are all entirely consistent with those <fO|[p=| 
r perovskite insulators such as BaTiOsE^I'EHl 
all of which possess rather typical lattice 
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dielectric constants of eph — 10-100. 

In summary, our first-principles calculations indicate 
that both CCTO and CdCTO possess similar intrinsic 
structural, vibrational, and dielectric properties. Our 
computations of the lattice contributions eph(O) to the 
static dielectric constant are in good order-of-magnitude 
agreement with the values measured experimentally in 
the far-infrared range for both CdCTO and CCTO, but 
are drastically smaller than the enormous values mea- 
sured at frequencies below the Debye cutoff range. The 
latter discrepancy reinforces the conclusion that some ex- 
trinsic mechanism is likely to be responsiblcJor the large 
dielectric constant present in both materials. u Several im- 
portant issues remain to be resolved by experiment, such 
as the connection between internal sample morphology 
and the large response, and the origin of the anomalous 
behavior of the oscillator strengths of the Cu-0 mode at 
low frequency observed in both materials. 
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